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ABSTRACT 



Type la supernova (SNIa) explosions synthesize a few tenths to several tenths of a solar mass, whose composition is the result of 
incomplete silicon burning that reaches peak temperatures of 4 GK to 5 GK. The elemental abundances are sensitive to the physical 
conditions in the explosion, making their measurement a promising clue to uncovering the properties of the progenitor star and of the 
explosion itself. Using a parameterized description of the thermodynamic history of matter undergoing incomplete silicon burning, 
we computed the final composition for a range of parameters wide enough to encompass current models of SNIa. Then, we searched 
for combinations of elemental abundances that trace the parameters values and are potentially measurable. For this purpose, we divide 
the present study into two epochs of SNIa, namely the optical epoch, from a few weeks to several months after the explosion, and the 
X-ray epoch, which refers to the time period in which the supernova remnant is young, starting one or two hundred years age and 
ending a thousand years after the event. During the optical epoch, the only SNIa property that can be extracted from the detection of 
incomplete silicon burning elements is the neutron excess of the progenitor white dwarf at thermal runaway, which can be determined 
through measuring the ratio of the abundance of manganese to that of titanium, chromium, or vanadium. Conversely, in the X-ray 
epoch, any abundance ratio built using a couple of elements from titanium, vanadium, chromium, or manganese may constrain the 
initial neutron excess. Furthermore, measuring the ratio of the abundances of vanadium to manganese in the X-ray might shed light 
on the timescale of the thermonuclear explosion. 

Key words, nuclear reactions, nucleosynthesis, abundances - supernovae: general - white dwarfs 



1. Introduction 



Type la supernovae (SNIa) are the brightest phenomena powered 
by nuclear reactions in the universe, releasing ~ 1.5 x 10 51 erg 
that are one-third invested in overcoming the binding energy of 
the progenitor white dwarf (WD), and the rest in kinetic energy 
of the debris. The nuclear energy release takes roughly one sec- 
ond, implying a mean power of ~ 1.5 x 10 44 W. These figures 
make SNIa prime targets for studies of nucleosynthesis in ex- 
treme conditions of density and temperature. 

While there is currently little doubt that the star explod - 
ing as a SNIa is a binary WD (e.g. Nugent et aD 1201 ll) . 
other details of the progenitor system are uncertain; to cite 
a few: the nature of the companion star (either double de- 



generate, DD, or single degenerate, SD; 


Patat etalJ 12007; 


Maoz & Mannuccil2008UBlondin et alJ2009t 


Simon et alJ200S : 


Li et alJ 201 It iBloom et alj|2012; Folev et al. 


20121 Maoz etal. 


2012: ISchaefer & Pagnotta 2012), the mass-accretion history 



120121). the geometry of the flame a t thermal runaway 



try ot tl 

dHoflich & Steidl200alZingale et al.ll2009l) . the mass of the WD 
at ignition (either a sub-Chandrasekhar W D, a Chandrasekhar- 
mass WD, or a super-Chandrasekhar one. | Mazzali et all [2007 ; 
Ivan Kerkwiik et ail 120101 : lYuan et all l2010l) . or the properties 
of the environment within a parsec of the WD. In the SD 
Chandrasekhar-mass WD scenario, the structure is supported by 
the pressure of degenerate electrons, which can be described by 
a few parameters. As a consequence, only a handful of the pro- 
genitor system properties have an impact on the thermonuclear 
explosion, namely the central density of the WD, the initial ge- 



ometry of the flame, and the chemi cal composition. Th e thermal 
content may also have an impact dCalder et al.ll2012h . Current 
knowledge of stellar evolution implies that the progenitor WD 
has to be made up mainly of carbon and oxygen, in propor- 
tions that depend on its main-sequence mass and metallicity. The 
metallicity has further influence on the WD structure through 
setting the electron mole number, Fe, which controls the num- 
ber of electrons available to support the WD. It is to be expected 
that nature makes SNIa from a range of all the above parameters. 

The workings of the thermonuclear bomb that powers a 
SNIa are not precisely known. Within the SD Chandrasekhar- 
mass WD scenario, the favored model involves an initially sub- 
sonic flame (deflagration) followed by a detonation. The release 
of thermonuclear energy in either burning-propagation mode is 
what controls the time schedule of the explosion. Thus, there is 
an intimate relationship between the final energy and mechanical 
structure of the supernova ejecta, on the one hand, and its chem- 
ical composition and profile, on the other hand. Furthermore, the 
chemical structure of the ejecta reflects the thermodynamic con- 
ditions in the explosion. Thus, there would be possible in prin- 
ciple to make an inverse analysis of the ejecta composition to 
uncover the details of the supernova explosion. 

From the observational point of view, obtaining the chemi- 
cal composition of the ejecta is not an easy task. Optical spec- 
tra depend in a complex way on the physical conditions in the 
ejecta days to weeks after the supernova explosion, which in 
turn are controlled by the energy deposition from radioactive 
nuclei, mainly 56 Ni and its decay product 56 Co. However, re- 
cently there have been advances that allow performing detailed 
tomographies of a few well-studied SNIa dStehle et al.l l2005a: 
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lAltavilla et aill2007t iMazzali et alJl200l iTanaka et alJl201 lh . It 
is to be expected that the set of SNIa for which a detailed chem- 
ical composition of the ejecta is available will increase steadily 
in the next several years. 

An alternative way to gain knowledge about the chemical 
struct ure of the ejecta is through studying yo ung supernova rem- 
nants (Ba denesl2010l) . lBadenes et al.l (l2008) studied the relation- 
ship between X-ray spectral features due to Mn and Cr in young 
SNIa remnants and the metallicity of the WD progenitor, which 
allowed them to conclude that the metallicity of the progeni- 
tor of Tycho's supernova is supersolar. The same technique has 
been since applied to other suspected SNIa remnants to deter- 
mine their progenitors' metallicity. 

The picture of SNIa ejecta drawn from observations has con- 
verged to a layered ch emical structure, alth ough with exceptions 
and peculiarities (e.g. iFolatelli et alJl2012h . In the center, there 
is a vo lume filled with re latively neutron-rich iron-group ele- 
ments (HoflichetaL||2004), thought to arise from matter heated 
to more than ~ 5.5 billion degrees which achieved nuclear sta- 
tistical equilibrium (NSE) and thereafter was neutronized effi- 
ciently by electron captures. This neutronization clears the cen- 
ter of the main radioactive product of SNIa, 56 Ni. Above the 
central core, there is a volume made of iron group elements that 
achieved NSE at a density low enough to experience a negligi- 
ble amount of neutronization. This volume is mainly composed 
of 56 Ni. Outwards from this region, there is a volume made of the 
products of incomplete silicon burning (hereafter, Si-b), cover- 
ing a wide range in atomic number, from silicon to iron. Finally, 
close to the surface of the ejecta, there is a tiny region made 
of unbumed fuel and material that has been processed no fur- 
ther than oxygen burning. The details of the layered structure 
vary from object to object, and not all the observations point to 
the same degree of chemical stratification. Some objects seem 
to have experienced a thorou gh mixing of the t wo outer layers, 
unburned and Si-b matter (Stehle et al. 2005b), although there 
is evidence that iron-group nuclei and the products of Si- b have 
clearly different physical histories (Rakowski et al. 2006). 

Although the masses of the different layers of SNIa ejecta 
are thought to vary from object to object, it can be said that 
NSE matter should amount from 0.2 to 0.8 M©, whereas 
intermediate-mass elements, mainly made as a result of Si-b, 
cover more or less the same range of masses (Mazzali et al. 
120071) . In principle, the study of both regimes therefore has the 
capability of providing interesting clues to the nature of the ex- 
plosion. However, whereas the abundances of NSE matter are 
influenced by the amount of electron captures achieved in the 
first phases of the explosion, the composition resulting from Si- 
b only depends on the initial composition of the supernova pro- 
genitor and on the dynamics of the explosion. 

The purpose of the present work is to elucidate the depen- 
dence of the chemical composition arising from Si-b on the 
progenitor properties (metallicity and chemical composition of 
the exploding WD) and on the explosion dynamics (expansion 
timescale and combustion wave velocity). We note that, although 
incomplete silicon burning has been studied in depth in quite 
a few works (e.g. Bodansk v et al] 119681: IWooslev et al] Il973t 
iHix & Thielemann||l996l. Il999t iMever et al]|1996lll998h . none 
of them have looked at the final chemical composition for trac- 
ers of the parameters governing this nucleosynthetic process. We 
use a model of the thermodynamic evolution of matter under- 
going Si-b, which we explain and justify in the next section. 
Thereafter, we focus on abundance ratios of Si-b products, which 
allows us to avoid the problem posed by the dependence of the 
absolute abundances on the total energy of the supernova. We 



split the results into two groups of abundance ratios: on the one 
hand, those that might create a detectable signature in the opti- 
cal display of the supernova, encompassing between a few days 
and a few months after the explosion, which we term the opti- 
cal phase or the optical epoch; and on the other hand, those that 
would apply to supernova remnants several hundred years after 
the event, which we term the X-ray phase or the X-ray epoch, 
because abundance measurements in this phase are usually ob- 
tained from X-ray spectra. We also refer to the first second of the 
explosion as the nucleosynthetic epoch, since it is then when the 
nuclear fluxes are able to modify the final chemical composition, 
apart from radioactive disintegrations. 

2. Numerical method 

We performed the nucleosynthetic calculations described in this 
paper by solving the nu clear kinetics equations w i th the code 
CRANK described in iBravo & Martmez-Pinedol d2012l) . For 
completeness, we summarize here its main features. CRANK 
integrates the temporal evolution of the abundances in a nuclear 
network for given thermal and structural (density) time profile, 
and for initial composition. The inputs to CRANK are the nu- 
clear data and the thermodynamic trajectories, as a function of 
time. The thermodynamic trajectories, in turn, can be the result 
of a hydrodynamic calculation of an explosive event, in which 
case CRANK provides the final chemical profile of the ejecta 
after computing the nucleosynthesis of each mass shell indepen- 
dently, without accounting for chemical diffusion. Alternatively, 
CRANK can compute the nucleosynthesis in a single homoge- 
neous region using an arbitrary thermodynamic trajectory. In the 
present work, we use CRANK in the latter way. 

T he nuclear evolutionary equations (see e.g. Rau scher et al .1 
|2002|) follow the time evolution of the molar fraction of each 
nucleus until the temperature falls below 10 8 K, after which 
time the chemical composition is no longer substantially mod- 
ified, with the exception of radioactive disintegrations. The in- 
tegration of the nuclear network follows the implicit iterative 
method of Wagoner with adaptive time steps. The iterative pro- 
cedure ends when the relative variation in the molar abundances 
is smaller than 10~ 6 , taking only those nuclei with molar abun- 
dance Y > 10~ 14 mol g _1 into account. The number of iterations 
is limited to 7, and if this limit is reached a new iteration pro- 
cedure with a reduced time step is started from the abundances 
at the last successful integration. In case the relative variation in 
the abundance of any nuclei with Y > 10~ 14 mol g _1 is greater 
than a prescribed tolerance, the iteration procedure is restarted 
from the last successful integration, but with a smaller time step. 
The sum of mass fractions is checked upon completion of each 
iteration procedure. If that sum differs from one more than 10 -8 , 
the iteration procedure needs to be restarted with a smaller time 
step. 

The network is adaptive, in the sense that its size (nuclei and 
reactions linking them) varies during the computation in order to 
improve its performance. A nucleus enters the calculation only 
if either it has an appreciable abundance (Y > 10~ 24 mol g -1 ) or 
it can be reached from any of the abundant nuclei by any one of 
the reactions included in the network. Light particles (neutrons, 
protons, and alphas) are the exception to the rule, because they 
are always included in the network. A reaction rate is included 
in the network only if the predicted change it induces on a molar 
abundance in the ne xt time step is larger t han a threshold, fixed 
at 10~ 20 mol g _1 (see Rausch er et al.l2002l for a similar method). 

The nuclear network consists of a maximum of 722 nuclei, 
from free nucleons up to 101 In, linked by three fusion reactions: 
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3a, 12 C+ 12 C, and 16 0+ 16 0, electron and positron captures, p~ 
and /3 + decays, and 12 strong interactions per each nucleus with 
Z > 6, involving neutrons, prot ons, alphas, and y. We tak e 
the neces s ary nu cl ear data fro m Salpeter & van Hornl (11 969): 
lltohetalj d 19791): iFuller et all (119821) : iMartmez-Pinedo et alj 
d2000l) : ICvburt et all d20ldl) ~ 



3. Modeling the physical conditions in incomplete 
silicon burning in SNIa 

Explosive nucleosynthesis studies can be divided into two broad 
categories, namely those that follow the thermodynamic tra- 
jectories computed from hydrodynamic models of the explod - 
ing object (e.g. iThielemann et al.l 119861: iMaeda et al.l [2010), 
and those that adopt a simplified parameterized description of 
the time dependence o f the relevan t thermodynamic variables 
(IWooslev et al.l 11973b iMever et al l 119981: Iffix & Thielemann! 
1 19991) . The former ones are more adequate for determining the 
chemical composition of the ejecta of a particular supernova (or 
whatever exploding object) model. The second ones are suitable 
for studies aiming at elucidating nucleosynthetic properties ex- 
pected to be shared by a general class of explosion models. Here, 
we chose this second approach because we are interested in just 
one nucleosynthetic process among all that are at work in SNIa. 

As is customary in parameterized studies of high- 
temperature explosive nucleosynthesis (op cit), we model the 
temporal evolution of the mass zone undergoing Si-b as a fast 
heating followed by adiabatic expansion and cooling. The tem- 
perature as a function of time is given by 

T _ | r «.o + ( r 9,peak _ r 9,«) '/frfse : < f < r rise 
9 1 r 9,peak ex P[-( f - T rise)/ T ] : r rise < f 

where t is time, Ty is the temperature in units of 10 9 K, and the 
rest of the quantities are model parameters, namely T r j se is the 
time to reach the maximal temperature, r 9 p ea k, starting from 
an initial value, T<)$, and r is the expansion timescale. In some 
studies, the expansion timescale has been taken as the free fall 
timescale, = 446/ y[p~\ however, we prefer to keep it as a 
free parameter in the present work. The initial temperature, T^o 
is irrelevant as long as it is low enough compared to the peak 
temperature. We fixed its value at — 1 in all the numerical 
experiments. 

The density is given by 



Pi = 



Pv,o + (p 7;P eak - P7,o) t /T rise : < t < r r 
P7,peak ex P[- 3 ( f - T rise)/ T l : T rise < f 



(2) 



where pi is the density in units of 10 7 g cm 3 , p 7 p ea i ( is the 
maximal density, and is the initial density, whose value is 
not influential in the present results and is fixed in the calcula- 
tions at pjfl = 0.3p 7 p ea j,. The factor three in the exponential 

in Eq. [2] accounts for the p cc T 3 dependence, appropriate to 
describing the isentropic evolution of a radiation-dominated gas 
dHix & Thielemannlll999l) . 

Although adiabatic expansion is a common assumption, mat- 
ter heated to temperatures high enough to achieve partial or total 
nuclear statistical equilibrium, T ^ 4 x 10 9 K, may go off the 
isentropic thermodynamic path because of the continuous ad- 
justment of the chemical composition as temperature declines, 
which releases nuclear energy. The degree to which Eqs. Q~]and[2] 



provide a faithful representation of the nucleosynthesis in SNIa 
layers undergoing Si-b can be tested, for any given SNIa model, 
by computing the chemical composition from these equations 
with parameters fitted to the thermodynamic trajectories pro- 
vided by the hydrodynamic code, then comparing it to the chem- 
ical composition obtained with the original thermodynamic tra- 
jectories. 

Figure [T] compares the yields obtained for the same ele- 
ments with the two thermodynamic time evolutions described in 
the previous paragrap h, f or the SNIa model DDTc describe d in 
iBadenes etakl d2005l) and lBravo & Martmez-Pinedol d2012l) and 
used hereafter as a reference model. As can be deduced from 
the figure, the errors in the abundances of elements from silicon 
to niquel related to the use of Eqs. [1] and [2] range from -40% 
to +70%. Most important, taking ratios of the yields of differ- 
ent elements, the errors can be made very tiny. For instance, the 
errors of the yields of the elements belonging to the first quasi- 
statistical equilibrium (QSE) group, i.e. those between silicon 
and scandium, are all very similar, implying that the ratio of the 
yields obtained with Eqs. Q] and [2] are a very good representation 
of the ratios belonging to the full thermodynamic trajectories 
provided by the supernova model. The errors of the elements 
belonging to the second QSE group, i.e. those between titanium 
and zinc, are much more heterogeneous. However, couples of el- 
ements can be identified in Fig.[T]whose yields are affected by a 
very similar error, hence for whom the yield ratios are predicted 
well by using Eqs. [T] and [2] An example is the couple formed by 
titanium and vanadium, at both one day and one hundred years 
after the explosion. Another example is the couple formed by 
chromium and manganese at one hundred years after the explo- 
sion. 

We also tried fitting the thermodynamic histories of super- 
nova layers with more complicated algebraic functions, but we 
obtained no improvements over the single exponential depen- 
dence on the adiabatic expansion timescale. 

3. 1 . Parameters ranges 

With the temporal evolution of temperature and density given 
by Eqs. Q] and [2] the full set of parameters that controls the final 
chemical composition of mass layers undergoing incomplete Si- 
b are the peak temperature T 9 pg^, the peak density p 7 p ea k, the 
expansion timescale t, the temperature rise timescale r r j se , the 
initial mass fraction of 12 C, X( l2 C), and the initial mass fraction 
of 22 Ne, X( 22 Ne). The last two quantities provide a convenient 
parameterization of the chemical composition of the WD at run- 
away. In this work, we simply assume that the only species ini- 
tially present in the WD are 12 C, 16 0, and 22 Ne. The initial mass 
fraction of oxygen is thus 

X( l6 0) = 1 - X( l2 C) - X( 22 Ne) , 
and the initial neutron excess is 
X( 22 Ne) 



(3) 



n = 



n 



(4) 



In turn, the initial neutron excess is related to the metallicity of 
the su pernova progenitor, Z, through rj » 0.101Z dTimmes et al.l 
120031) . which leads to, Z * X( 22 Ne)/l.llQ 

In the following, we discuss the relevant ranges of the pa- 
rameters to explore. 

1 It should be recalled that the relationship between the neutron ex- 
cess of the progenitor WD at runaway and the metallicity of the progen- 
itor star at zero-age main-sequence can be modified owing to electron 
captures during the hundreds of years of the carbon-simmering phase 
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Fig. 1. (Color online) Error in the elemental yields obtained with 
Eqs. \T\ and [2] relative to those for the same element obtained 
with the thermodynamic trajectories provided by a hydrody- 
namic code for a SNIa model (see text for further details). Large 
open (red) circles give the errors belonging to a time of one day 
after the explosion, when the abundances of many iron-group el- 
ements are still dominated by radioactive isotopes. Small filled 
(blue) circles give the errors belonging to a time one hundred 
years after the explosion. Points belonging to some interesting 
couples of elements are highlighted by a small ellipse around 
them. For this comparison, we have only considered the layers 
with peak temperatures between 4.3 and 5.2 GK. 

3.1.1. Peak temperature 

In ou r reference model, DDTc, as well as in most SNIa models 
(e.g.lWooslev et al.ll973l:rHix & Thielemann! 1999HMaeda et alj 
2010), incomplete Si-b is ac hieved at densities slightly in ex cess 
of 10 7 g ctrT 3 (see Fig. 1 in lBravo & Martfnez-Pinedol2012l) . At 
these densities, to achieve the QSE of the iron group and of the 
silicon group, the peak temperature has to exceed ~ 4 x 10 9 K. 
On the other hand, to avoid reaching full NS E the tempera- 
ture s hould remain well below ~ 5.5 x 10 9 K dWooslev et alj 
Il973h . As we demonstrate later, the most promising targets of 
the present study are the lightest elements of the iron-group, 
from titanium to manganese. A s can be seen in Table 2 of 
Bravo & Martinez-Pined^ (1201 2l) . the bulk of these elements 
originates in SNIa in layers attaining temperatures in the range 
from 4.2 x 10 9 K to 5.2 x 10 9 K, with the exception of man- 
ganese, which has a sizeable contribution from layers reaching 
NSE. For the purposes of the present study, we define the range 
of peak temperatures able to provide incomplete Si-b such as 
that between 4.3 - 5.2 x 10 9 K. 

Figure [2] shows the evolution of the nuclear species com- 
puted with r 9peak = 4.8, p 7p eak = 2 - 2 ' T = 029 s > T rise = 
10~ 9 s, X( l2 C) = 0.5, and X( 22 Ne) = 0.01. For this set of param- 
eters, and in spite of the high temperature attained, there is no ap- 

that precedes thermal runaway (Piro & Bildsten 2008; Badenes et al. 
2008). 



Fig. 2. Evolution of the chemical composition for a reference 
set of parameters, namely - 4.8 10 9 K, T r j se = 10~ 9 s, 

r = 0.29 s, X( 12 C) = 0.5, X( 22 Ne) = 0.01. The dashed line is the 
temperature in 10 9 K. 



preciable change in the chemical composition before t = r r j se . 
Carbon is exhausted in a time of a few tens of ns and oxygen 
burns at less than 10 fis, whereas the isotopes belonging to the 
silicon QSE-group achieve stable proportions after ~ 0.1 ms. It 
can be seen that the species of the iron QSE-group steadily in- 
crease in abundance after ~ 1 ms and until the expansion cools 
matter below ~ 4 x 1 9 K. After this point, two different trends in 
the isotopes belonging to the iron QSE-group can be discerned. 
On the one hand, 56 Ni is the only species that continues growing 
in abundance. On the other hand, the mass fractions of the rest 
of the most abundant components of the iron QSE-group, rich in 
neutrons, decrease owing to recombination with excess free pro- 
tons and alphas, adapting their abundances to equilibrium at the 
decreasing temperatures, w hile they remain above T > 3 x 10 9 K 
dHix & Thielem ann 1999). This recombination causes a drop in 
the abundance of 4 He after ~ 0.01 s, which is clearly visible in 
the plot. 

It is notable that the final mass fractions of 55 Co ( 55 Mn after 
a few years of radioactive disintegration) and 52 Fe (grandparent 
of 52 Cr) are among the largest from the iron QSE-group. But for 
54 Fe and 56 Ni, they are the most abundant species with 77 + 
and with 77 = 0, respectively. 

3.1.2. Distribution function of peak temperatures. 

In any model of SNIa, matter attains a wide range of peak tem- 
peratures. As the burning front moves through layers of decreas- 
ing density, its peak temperature drops. Matter experiencing Si-b 
is not subject to just a single peak temperature but to a sequence 
of them. Figure [3] shows the distribution of peak temperatures 
in a suite of SNIa models of different flavors. The distribution 
function of peak temperatures in these models can be described 
by either an approximately power-law dependence on T 9 pg^, 
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which we fit with dm/dlpg^ oc T^^, or as independent of 

T 9 pg^- Aside from the models shown in these figures, model 
O-DDT in lMaeda et alj ( 1201 fj) follows a dependence on T 9 
that is similar to the proposed power law, at least qualitatively, 
whereas those of models W7 and C-DEF (op cit) are described 
better as independent of T 9 peak- 
In the following, we use these two representative distribution 
functions of the peak temperatures to estimate the mass of the 
different elements created by given combinations of the param- 
eters of Si-b evolution, and we analyze the impact of the peak- 
temperature distribution function on the results. If this impact 
were sizeable, observations might provide insight into the dis- 
tribution function of peak temperatures in SNIa. If, on the other 
hand, its impact is negligible, the results will be independent of 
the assumed distribution function. 



3.1.3. Peak density 

Peak density in Si-b depends on the peak temperature, both 
quantities being correlated with entropy. iMever et aD (ll998l) in- 
troduced a quantity proportional to the entropy per nucleon in 
photons, 



10 



ee 3.4 x 10" 3 r 9 3 /p 7 



(5) 



Given the thermodynamic history described by Eqs. Q]and[2] <p 
remains constant. 

For model DDTc, in the region undergoing incomplete Si-b 
the peak density, p 7 pg^ can be approximated as a function of 
the peak temperature by the relationship 



P7,peak = ex p( r 9,peak ~ 4 ) 



(6) 



Using this equation, <p ranges fro m 0.15 to 0.2 , which are much 
lower values than investigated in Meveret al. (1998). In SNIa, 
the peak temperature and density are the result of the thermo- 
dynamic conditions in matter that is either shocked in a deto- 
nation wave (delayed- detonation mode ls) or burnt subsonically 
(deflagration models). Khokhlov ( 1988) provides the peak tem- 
perature and density for Chapman- Jouguet deflagration and det- 
onation waves, which differ from Eq. |6]by at most 22% in the 
temperature and density ranges of Si-b. 

To reduce the number of free parameters, we use Eq. |6]in 
the calculations reported in the present work. It is therefore nec- 
essary to evaluate the degree of uncertainty derived from this 
relationship. This is done in Fig. [4] where we represent the ra- 
tio from manganese to chromium at the X-ray epoch, for peak 
densities given by either Eq. [6] a factor of two higher or a fac- 
tor of two lower. It is clearly seen that the precise value of the 
peak density is irrelevant to the mass ratio of these elements. The 
points belonging to the different peak densities used are nearly 
indistinguishable in the plot, except at the highest neutron ex- 
cess and the lowest temperature. These last points belong to the 
highest ent ropies studied in th e present work. It is interesting 
to note that M eyer et al.l d!998l) explored even highest entropies 
and neutron excesses. 



3.1 .4. Rise time 

Most, if not all, nucleosynthetic studies of the explosive incom- 
plete Si-b process compute the nuclear flows st arting from the 
peak t emperature and the initial fuel composition. IWooslev et al] 
(I1973I) argue that this is justified when the rise time, i.e. the time 
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Fig. 4. (Color online) Final mass fraction ratio Mn/Cr, at a time 
of 100 yrs after the explosion, as a function of peak temperature 
in incomplete Si-b for different maximal densities (point type) 
and initial neutron excesses (color online). The maximal densi- 
ties are given by either Eq. [6] (filled circles), or by a factor of 
two higher (triangles), or lower (crosses) than Eq. |6] The initial 
neutron excesses are rj = 2.27 10~ 4 (red), 9.09 10~ 4 (blue), and 
6.82 10~ 3 (black), increasing from bottom to top. 



taken to dump the energy needed to raise its temperature to its 
peak value into the mass zone, is much shorter than the expan- 
sion timescale, r r j se <k t. In SNIa models, the energy dumped 
into a zone until it reaches its peak temperature can either be due 
to a compressional heating related to the passage of the shock 
front at the leading edge of a detonation wave or it is a combi- 
nation of thermal diffusion and energy released by nuclear reac- 
tions in a deflagration wave. The rise time associated with each 
one of these burning waves is very different. 

Figure shows a comparison of both rise times as a func- 
tion of density, together with the hydrodynamic timescale. The 
thermal structure of a laminar flame (deflagration) is determined 
by heat diffusion from hot ashes to cool fuel at low temperature 
and by the release of heat by nuclear reactions at high tempera- 
ture. The transition between both regimes takes place at a criti- 
cal temperature, T cr j f , that depends on density. The deflagrative 
timescale in Fig.|5]nas been estimated as the nuclear timescale at 
•^crit' F° r simplicity i t> om T CI ^ and the nuclear timescale were 
computed, accounting only for the dominant reaction rate, i.e. 
the fusion of two 12 C nuclei. 

Triggering of ignition in a supersonic combustion wave (det- 
onation) is due to an initial temperature jump in a length scale 
of a few mean free paths, in a shock front. Later on, the heat 
released by thermonuclear reactions determines the thermal pro- 
file. Thus, the detonation timescale can be estimated as the nu- 
clear timescale at the shock temperature, T^j™^. For the cal- 
culations shown in Fig. [5] the shock temperature was calculated 
by solving the Hugoniot jump conditio ns using the veloc ities of 
Chapman- Jouguet detonations given in Khokhlov! lfl988h . 
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Fig. 3. (Color online) Distribution of maximal temperatures in the incomplet e Si-burning region obtained from a suite of SNIa 
models. Left: One-dimensional delayed detonati on model s DDT c (red crosses. iBadenes et alj|2005l) and DDTf (red solid squares), 
a sub-Chandrasekhar model ( blue stars) fr om Bravo et al.l (1201 ll fourth row in their Table 2), and three-dimensional delayed det- 
onation model TURB7 from iBravo et al.1 (12009 , red empty squares,). The solid line is a fit to the data given by, dm/dTpg^ oc 

T n , . Right: Three-dimensional delayed detonation models DDT3DA (red empty squares) and DDT3DB (red crosses) from 

Brav o & Garcfa-Senzl d2008). gravitationally confined deto nation mode l GC D1 (magenta solid triangles), pulsating reverse deto- 
nation model PRE ) 18 (m agenta empty triangles), both from lBravo et al.l d2009l) . and sub-Chandrasekhar model C (blue stars) from 
iGarcfa-Senz et al.l (l999). The distribution of maximal temperatures in these models is approximately constant, independent of the 
peak temperature, unlike the models shown in the left panel. 



The detonation rise time is the shortest one by orders of 
magnitude, whereas the deflagration rise time lies between the 
detonation rise time and the hydrodynamic timescale. An alter- 
native way of estimating the timescale of a deflagration is as 
the ratio between the flame width and velocity. U sing the flame 
properties reported in Timmes & Wooslev (1992), the deflagra- 
tion timescale goes into the range 9 x 10~ 3 - 10~ 4 s for densities 
between 10 7 and 5 x 10 7 g crrT 3 . Although these timescales are 
about one order of magnitude longer than the ones in Fig. [5] the 
qualitative conclusions remain the same. 

In view of the large difference between both rise times and 
the hydrodynamic timescale, it is to be expected that the pre- 
cise value of T r j se will be unimportant for the nucleosynthesis. 
Notwithstanding this expectation, we have explored the nucle- 
osynthesis resulting from different rise times ranging from 10 -9 s 
to 0.2 s, in search for possible, although improbable, tracers 
of the type of burning wave responsible for incomplete Si-b in 
SNIa. 

Figure [6] shows the evolution of nuclear species for condi- 
tions similar to those in Fig. [2] apart from the rise time. In the 
left hand panel of this figure, the rise time has been set to 10~ 5 s 
to illustrate the case of a deflagration wave, whereas Fig. |2]is 
more representative of a detonation wave. In spite of some dif- 
ferences in the evolution of the chemical abundances, as for in- 
stance in the extent of burning before Tg pg^, the final mass 
fr actions in both figures are identical, in line with the predictions 
of lWooslev et al.l(ll973l) . Moreover, the same conclusion applies 



to a calculation with a rise time, r r j se = 0.1 s, comparable to the 
expansion timescale (right panel in Fig. 

3.1.5. Expansion timescale 

In SNIa, the expansion timescale is bounded to be longer than 
the WD sound-crossing time, r > ? sounc j = ^WD/ v sound ~ 

0. 1 - 0.2 s, which is similar to the free-fall timescale at the den- 
sities of interest to incomplete Si-b. In the present work, we ex- 
plore expansion timescales in the range 0.1 - 0.9 s. 

The left hand panel of Fig. [7] shows the evolution of the 
chemical composition for an expansion timescale t = 0.72 s, 

1. e. slower than the one in Fig. [2] r = 0.29 s, while the rest of 
parameters remain unchanged. As could be expected, the differ- 
ences only appear in the final phase of the evolution, as can be 
seen in the temperature curve and in the leveling out of the abun- 
dances when all the species leave QSE. Interestingly, a longer 
expansion timescale favors synthesis of a larger quantity of 56 Ni 
at the expense of 28 Si and 32 S, but the rest of the species with 
X > 0.001 is hardly affected at all. 

3.1.6. Initial carbon abundance 

Domfngue z et al.l (l200lh accurately computed the evolution of 
stars in the range 1.5-7 M0 from main sequence to the for- 
mation of a WD, starting from different metallicities. According 
to their results, the carbon mass fraction in the WD lies in the 
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Fig. 6. Evolution of the chemical composition for a rise time to the maximal temperature different from that in the reference set 
of parameters (Fig. [2}. Left: Slower rise to the peak temperature, T^ se = 1CT 5 s. Right: Rise time comparable to the expansion 
timescale, r r j se =0.1 s. 
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Fig. 7. Evolution of the chemical composition for different expansion timescales and initial carbon abundances than in the reference 
set of parameters (Fig. |2j. Left: Slower expansion timescale, t = 0.72 s. Right: Smaller carbon abundance, X( l2 C) = 0.3. 

range 0.37 - 0.55. We therefore explore a slightly wider range elements. Once oxygen is exhausted, the evolution in both fig- 
fromX( 12 C) =0.3toX( 12 C) = 0.7. ' ' ' ures is identical. 



The right hand panel of Fig. Q shows the evolution of the 3.1.7. Progenitor metallicity 
chemical composition for initial abundances X( l2 C) = 0.3 and 

X( 16 0) = 0.69. Comparison with Fig.EJ in which X( 12 C) = 0.5 Type la supernovae are detected in host galaxi es with a wide 

and X( 16 0) = 0.49, reveals that the initial carbon abundance has range of metallicities (e.g. Z « 0.1 - 2.5Zq in ID' Andrea et al.l 

no influence on the final nucleosynthesis of the most abundant 201 1), with apparent preference for hosts with Z > Zq. The low 
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Density (10 7 g cm" 3 ) 



Fig. 5. Comparison of relevant timescales as a function of the 
initial density. Solid line: expansion timescale, computed as the 
free-fall or hydrodynamic timescale, 446/ y[p s; dotted line: rise 
time in a deflagration wave, computed as the nuclear timescale at 
the critical temperature at which the nuclear timescale matches 
the heat diffusion timescale; dashed line: rise time in a detona- 
tion wave, computed as the nuclear timescale at the temperature 
of shocked matter. 



end of the metallicity range migh t reaches Z » 0.04Z (T) accord- 
ing to the preliminary results in Riga ult et al.l (1201 ll) . Here, we 
present results of incomplete Si-b for initial 22 Ne mass fractions 
between 2.5 x 1(T 4 and 7.5 x 1(T 2 i.e. Z « 0.02 - 5Zfp for a so- 
lar me tallicity Zp * 0.012-0.014 dLoddersl2003HAsplund etall 
120061) . 

Figure [8] shows the evolution of the chemical composition 
starting from a neutron excess 7.5 times longer than in Fig. [2] 
Comparing these two figures, the impact of the initial neutron 
excess on the final abundances is evident. In Fig. [8] the initial 
amount of 22 Ne disappears as soon as carbon is consumed. The 
role of 22 Ne as the main neutron-excess holder is later played by 
34 S in oxygen burning, and by 54 Fe and 55 Co in silicon burning 
and during freeze-out. 

3.2. Summary of model parameters 

To summarize, four free parameters remain for which we try 
to find tracers from the nucleosynthesis of incomplete Si-b, 
namely the temperature rise time, the expansion timescale, the 
initial carbon abundance, and the progenitor metallicity, Z m 
X( 22 Ne)/l.l 1. The initial values of temperature and density, Tg^ 
and pi o, are fixed, and the peak density is linked to the peak tem- 
perature through Eq. |6] Finally, the peak temperature is another 
parameter that is varied in the different nucleosynthetic calcula- 
tions, but such results are integrated using the distribution func- 
tion of peak temperatures. To be precise, in the following when 
we present the results for a given set of parameters [T r j se , r, 



Fig. 8. Evolution of the chemical composition for an initial 
neutron excess higher than in the reference set of parameters 
(Fig.H,X( 22 Ne) =0.075. 

Table 1. Sets of parameters of Si-b shown in Figs. [9] QT] and 

EHBJ'. 



Parameter Values 





11,0 


T 


(0.1,0.2, 0.3,0.4, 0.5, 0.6)/ln(2) s 


X( 22 Ne) 


2.5 x 10~ 4 , 2.5 x 10~ 3 , 0.01, 0.025, 0.075 


r rise 

X( r2 C) 


10- 9 , 10~ 7 , 10~ 5 , 10~ 3 , 0.01, 0.05, 0.1, 0.2 s 


0.3, 0.4, 0.5, 0.6, 0.7 



Notes. (n) All the combinations of the parameters values are represented 
in these figures except when the contrary is stated in the corresponding 
captions. {b) Exponent of Tp eg ^ in Eq. [JJ The integral in Eq. [7] was 
approximated by summing the abundances obtained at different tem- 
peratures, Tg = 4.3 - 5.2 with interval 0.05 GK, where each point was 

weighted by T" , . 
° J peak 

X( l2 C), and Z], the mass fraction of each element is given by 

x =/^ r peak)( dm / dr peak) d7 peak 

<*/X(7peak)7;eak d7 peak> V) 
with dm/dT^ ea ^ oc 7^^, and either n- 1 1 or n = (Fig. [3}. 

4. Tracers of the progenitor metallicity: X-ray epoch 

Since isotopic abundances cannot be measured in supernova 
ejecta, it is convenient to know what isotopes are the main con- 
tributors to the abundances of elements made in incomplete Si-b, 
in the X-ray as well as in the optical epochs. This information is 
given in Table [2] When focusing on the X-ray phase, all the ra- 
dioactivities that synthesize sizeable quantities of elements from 
scandium to iron cease a long time before forming a supernova 
remnant, so their elemental abundances remain constant at this 
epoch. 
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The only pairs of elements whose mass ratio in the X-ray 
epoch correlates strongly with progenitor metallicity are man- 
ganese vs chromium, manganese vs titanium, and vanadium 
vs titanium. We discuss them next, starting with the ratio of 
manganese to chro mium, which was originally proposed by 
Bade nes et alj d2008l) as a metallicity tracer. 



4.1. Ratio of manganese to chromium 

Figure [9] shows the ratio of the abundances, either total masses or 
mass fractions, of manganese to chromium (Mn/Cr) in the X-ray 
epoch, averaged by applying Eq. [7] as a function of the initial 
neutron excess. The insensitivity of the Mn/Cr ratio to the tem- 
perature rise time and to the initial carbon abundance is remark- 
able, as well as to the distribution function of peak temperatures. 
At any given neutron excess, the final Mn/Cr ratio varies by less 
than a factor of two for the whole range of the rest of parameters, 
whereas the Mn/Cr ratio as a function of 77 varies as much as a 
factor of 80. Furthermore, manganese and chromium are one of 
the element pairs whose ratio is determined better with the sim- 
plified thermodynamic evolution provided by Eqs. Q] and [2] (see 
Fig. [TJ. Thus, the Mn/Cr ratio stands out as a robust tracer of the 
initial neutron excess at WD runaway. The Mn/Cr ratio can be 
fit as a function of 77 by a power lavv0 (Fig. |9), 



M(Mn)/M(Cr) = 49.177' 



0.676 



(8) 



The results presente d here are complementary to the work of 
Bade nes et alJ d2008l) . where the Mn/Cr ratio is computed for a 
series of SNIa models. We find a similar power-law exponent as 
the one reported in that work, 0.676 vs 0.65, and t he new power- 
law fit stays within a factor of two from the one in Ba denes et al.l 
(2008). 

The robustness of the Mn/Cr ratio as a metallicity tracer is 
enhanced by its relative insensitivity to the peak temperature 
(Fig. [4]). As a result, the Mn/Cr ratio is insensitive to the degree 
of mixing the nucleosynthetic yields for the supernova layers un- 
dergoing incomplete Si-b. 

One can wonder if the Mn/Cr ratio can be used as a metallic- 
ity tracer for core-collapse supernovae. The answer is no, as long 
as layers experiencing Si-b in core collapse supernovae evolve 
through a much higher isoentropic than in SNIa. We show in 
Fig. [9] the results of the Mn/Cr mass ratio as a function of 77 
when the parameter <p defined in Eq.[5] takes the va lue <p _ 6.8, 
which is appropriate for core collapse supernovae dMever et al.l 
119961) . while t varies in the range 0.29 - 0.87 s, and the range of 
X( I2 C) is 0.3 - 0.7. It can be deduced from this figure that the 
tight relationship between the abundances of Mn and Cr and the 
metallicity is a result of a low entropy evolution that cannot be 
extrapolated to core-collapse conditions. 

4.1.1. Understanding the link between neutron excess and 
the final ratio of manganese to chromium 

It is instructive to try to understand the origin of the tight depen- 
dence of the Mn/Cr ratio in the X-ray epoch on the initial neu- 
tron excess. For this purpose, it is necessary to tell the story of 
the synthesis of their parent species, 55 Co and 52 Fe (see Table[2]i. 

At the peak temperatures studied here, both 55 Co and 52 Fe 
attain abundances in statistical equilibrium with the iron QSE- 
group. As temperature drops, both nuclei evolve to maintain 



2 The statistical uncertainty in the exponent of Eqs. [8] and [20] is 
±0.001, that of Eq.TOlis ±0.003, and that of Eqs. [22j [23] and l24l is 
±0.01. 




0.0001 



0.001 



0.01 



Fig. 9. (Color online) Averaged final mass fraction ratio Mn/Cr 
as a function of the initial neutron excess. The average mass 
fractions are drawn either from the fit to the distribution of maxi- 
mum temperatures given in the left panel of Fig.[3](red circles) or 
from a uniform distribution function representative of the models 
shown in the right panel of the same figure (blue crosses). In both 
cases, the results are shown for different r r j se in the range 10 9 - 
0.2 s, t in the range 0.10-0.90 s, and X( 12 C) in the range 0.3- 
0.7. The solid line is a fit given by, M(Mn)/M(Cr) = 49 Arf- 616 . 
The dashed line shows the relationship expected from an analytic 
model, Eq.[l7] Green stars give the mass fraction ratio Mn/Cr for 
high-entropy conditions (see text for details). 



this equilibrium, which is broken when T < 3 x 10 9 K 
dHix & Thielemann! 19991) . The left hand panel of Fig.fTOlshows 
the evolution of the abundances of 55 Co and 52 Fe during expan- 
sion and cooling for six different sets of parameters, given in 
Table [3] It can be seen that the abundance of 55 Co levels off 
at T - 2.5 x 10 9 K, whereas that of 52 Fe stabilizes even ear- 
lier, at T - 3.5 x 10 9 K. The highest differences between the 
curves belonging to 55 Co are found for the extremes metallicities 
tested: the bottom solid curve, X( 22 Ne) = 2.5 x 10~ 4 , and the top 
solid curve, X( 22 Ne) = 0.075. The abundance of 55 Co remains 
in equilibrium with theirs neighbors in the iron QSE-group until 
T - 2 x 10 9 K, as can be seen in the right hand panel of Fig. [10] 
In this figure, the molar fluxes of the main reactions linking 55 Co 
with other members of the iron QSE-group are plotted, where the 
most important are 54 Fe -Hp ±^ 55 Co +y and 56 Ni +y ±5 55 Co -Hp. 
The flux rate of the proton capture on 55 Co begins to depart from 
the rate of the (y, p) reaction on 56 Ni just after T 2.8 x 10 9 K. 
However, the net rate of destruction of 55 Co is always at least an 
order of magnitude lower than the rate of each of the reactions 



Fe + p 



Co + y, thus ensuring that equilibrium is main- 



tained in the temperature range shown in Fig. [TO] Assuming that 
the abundances within the iron QSE-group remain in equilib- 
rium allows us to formulate a simple model for the ratio of the 
final (before radioactive disintegrations) mass fractions of 55 Co 
to 52 Fe, X( 55 Co)/X( 52 Fe). 
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Table 2. Parent species of Si-b elements at the optical and X-ray epochs. 



Element 


Optical epoch" 


Variation 6 


X-ray epoch 1 ' 


Sc 


4 Ti (0.18 s) 




43 Ti (0.18 s) 


1 1 


cr (zi. o nj — > v (lo.u a) 


t 
1 


cr (zi.o nj — > v (lo. u 


V 


48 Cr (21.6 h) -> 48 V (16.0 d) 


1 


51 Mn(0.05 s)-> 51 Cr (27.7 d) 




49 Cr (0.04 s) -> 49 V (330 d) 








51 Mn (0.05 s) -» 51 Cr (27.7 d) 


T 




Cr 


52 Fe (8.3 h) -> 52 Mn (5.6 d) 


T 


52 Fe(8.3 h)^ 52 Mn (5.6 d) 


Mn 


52 Fe (8.3 h) -> 52 Mn (5.6 d) 


1 


55 Co (17.5 h) -> 55 Fe (2.7 y) 




53 Fe (0.009 s) -> "Mn (3.7 Myr) 






Fe 


54 Fe 




56 Ni (6.1 d) -> 56 Co (77 d) 



Notes. (fl) Main contributor to the abundance of the element at the optical epoch (in parenthesis, if it is an unstable isotope, its radioactive lifetime). 
<b) This column indicates if the contribution of a radioactive chain to the element abundance in the optical phase increases with time (f), decreases 
(|), or remains about constant (~). (r) Main contributor to the abundance of the element at the X-ray epoch (in parenthesis, if it is an unstable 
isotope, its radioactive lifetime). 



T | y i ^ 




Temperature (10 9 K) Temperature (10 9 K) 

Fig. 10. Nuclear evolution as a function of temperature during expansion and cooling. Left: Evolution of the mass fractions of 55 Co 
(solid lines), 52 Fe (dashed lines), and a (dot-dashed lines), for six different sets of parameters, as given in Table [3] The dotted lines 
show the quantity /, defined in Eq. [16] Right: Evolution of the molar fluxes of the main reactions contributing to the abundance 
of 55 Co for the same set of parameters as in Fig. |2] Reactions building 55 Co are drawn as solid lines, reactions destroying it are 
drawn as dashed lines, and the absolute value of the net rate of d5 / ( 55 Co)/df is drawn as a dotted line. Upper curves are for the 
reactions 54 Fe + p ±^ 55 Co + y, in which both line types are nearly indistinguishable. The middle curves are for the reactions 
56 Ni + y ±+ 55 Co + p. Finally, the bottom curves are for the reactions 55 Fe + p ±5 55 Co + n. 



Table 3. Sets of parameters of Si-b represented in the left panel 
ofFig.QI)] 



^9,peak 


r 


AT( 22 Ne) 


T rise 


X( ll C) 




(s) 




(s) 




4.8 


0.29 


0.010 


io- y 


0.5 


4.3 


0.29 


0.010 


10" 9 


0.5 


5.2 


0.29 


0.010 


10- 9 


0.5 


4.8 


0.72 


0.010 


10" 9 


0.5 


4.8 


0.29 


0.075 


10" 9 


0.5 


4.8 


0.29 


2.5 x 10" 4 


10-" 


0.5 



We begin by assuming that the initial neutron excess in the 
matter subject to incomplete Si-b is stored in the species 54 Fe 



and 55 Co at the end of the nucleosynthetic epoch, which can be 
justified in light of Figs.|2]and Figs.|6]-[8] Then, using Eq.|4j one 
can write 



2T22.0 - 2T 5 4 + T55 , 



(9) 



where ^22,0 = X( 22 Ne)/22 is the initial molar fraction of 22 Ne 
and for simplicity, we write the molar fraction of species A Z as 
Ya, in mol/g. Another relationship between the abundances of 
54 Fe and 55 Co can be derived from the assumption of equilib- 
rium of the forward and backward (p, y) reactions linking both 
nuclei, and the ones linking 55 Co with 56 Ni, 



f 55 



T54T56 



r 54py/l56yp 
f 55py^55yp 



1/2 



T54T56 



C( 54 Fe)C( 56 Ni) 



1/2 



,(10) 
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where r 54 p y is the rate of radiative proton capture on 54 Fe, A 56y p 4.2. Ratio of vanadium to titanium 



is the rate of photodisintegration of 56 Ni through the proton 
channel, and so on. The last equality is a result of the assump- 
tion of global equilibrium within the i ron QSE-group, where th e 
quantity C( A Z) is defined by (see, e.g. lHix & Thielemann! 1 996), 



B( A Z) 



kT 



(11) 



where G( A Z) is the partition function, B( Z A) the nuclear bind- 
ing energy, Np^ Avogadro's number, k Boltzmann's constant, and 



: 5.9417 x 10 33 r 9 3/2 . 



Solving Eq. QJJ for F54 and substituting in Eq. [9] gives the 
following relationship between the final abundance of 55 Co and 
the initial abundance of 22 Ne, 



2t 2 

755 = T 



1 + 1 



« 2 



where 21 has been defined as 



91 = 8F 5 ^ 2 exp 



12.19 



(12) 



(13) 



T 5 2 



F 56 C( S2 Fe)C(a) 



Y a C( 56 Ni) ' 
From this equation and Eq.[T2]one obtains 



X( 55 Co) 
X( 52 Fe) 

where 



/ 



1 + 



16F 2 2,( 
2I 2 



1 



(14) 



(15) 



/ ee 2.4 x 10 



. w pY a /68.4 



is shown in Fig.QJjto only vary slightly for 3.8 > Tg S 2.2 and to 
be nearly independent of the parameters of the Si-b nucleosyn- 
thetic calculation. Substituting approximate values for / 0.07 
and 21 0.011, one finally obtains the desired relationship be- 
tween the abundances of 55 Co and 52 Fe, and the neutron excess, 



X( 55 Co) 
X( 52 Fe) 



^ 0.07 Ul + 6.61 x I0 4 // - ! 



(17) 



This equation nicely fits the results of the nucleosynthesis cal- 
culations, as can be seen in Fig. [9] We stress that the approxi- 
mations needed to derive Eq. [17] imply that it is not necessarily 
better than the simpler power-law fit given by Eq. [8] which is 
shown in the same figure. However, its derivation allows us to 
discern where the dependence of M(Mn)/M(Cr) on the progen- 
itor metallicity comes from, and what its conditions of validity 
are. It also shows that any effect of the rest of parameters is of 
second order in comparison with that related to a variation in Z. 
Finally, it also allows us to understand that the power-law index 
of 0.676 is the result of the comparable magnitude of X( 22 Ne) 
and 91. 



The ratio of abundances of vanadium to titanium, V/Ti, provides 
an alternative way of measuring the initial neutron excess of the 
progenitor of a SNIa in the X-ray phase. The left hand panel of 
Fig. [TT] shows that V/Ti traces nicely the initial neutron excess 
apart from 77 55 10~ 4 , for which there is a substantial dispersion 
of this ratio. The best power-law fit is 



M(V)/M(Ti) = I4.877' 



0.726 



(18) 



Notice the similarity between the exponent of the power laws in 
Eqs. [8] and [18] According to Fig.Q] the ratio V/Ti obtained with 
the analytic model of Si-b given by Eqs. Q] and [2] matches the 
ratio in hydrodynamical simulations of SNIa very well. 

In the X-ray epoch, vanadium and titanium come from 5 1 Mn 
and 48 Cr, respectively (see Table |2}, whose abundances in the 
nucleosynthesis epoch are linked by a captures to those of 55 Co 
and 52 Fe. Using the equilibrium relationships, as in section lT. 1.11 
the abundance ratio of 51 Mn to 48 Cr can be shown to be propor- 
tional to the ratio of 55 Co to 52 Fe, with 



I( 51 Mn) 

Equation [T2l sets the functional dependence of 5^5 on X( 22 Ne) Jf( 48 Cr) 
(or 77): for X( 22 Ne) <K 2l 2 it leads to T 55 oc X( 22 Ne), whereas for 
X( 22 Ne) » 2I 2 the result is Y 55 oc X( 22 Ne) 1/2 . Since 91 ~ 0.011 
when global QSE breaks down (Tg 2.9), it turns out that we 
are in an intermediate regime, in agreement with the power-law 
exponent in Eq. [8] 

In QSE, the molar abundance of 52 Fe can be written as 



~ 76c::pf 3.146\ *( 55 Co) 
-0.76 exp ^ jr j x( s2 Fe) 



(0.15-0.30) 



X( 52 Fe) ' 



(19) 



i.e., the abundance ratio V/Ti is expected to range about 15-30% 
of the ratio Mn/Cr. Comparison of the left hand panel of Figs. [TT] 
to[9]confirms this approximate relationship. 



4.3. Ratio of manganese to titanium 

The ratio of abundances of manganese to titanium, Mn/Ti, pro- 
vides an alternative way of measuring the initial neutron excess 
in the X-ray phase. The right hand panel of Fig. [TTI shows that 
Mn/Ti nicely traces the initial neutron excess in the whole 77 
range explored here. The best power-law fit is 



(16) M(Mn)/M(Ti) = 1170t/' 



0.677 



(20) 



which is again similar to the exponents of the power laws in 
Eqs. [8] and [18] According to Fig. [TJ the ratio Mn/Ti obtained 
with the analytic model of Si-b given by Eqs. [TJ and [2] shows a 
relative error of approximately 20% relative to the ratio in hy- 
drodynamical simulations of SNIa. However, this error is much 
lower than the range of variation in M(Mn)/M(Ti) as a function 
of 77. 

In the X-ray epoch, manganese and titanium come from 55 Co 
and 48 Cr, respectively (see Table Using the equilibrium re- 
lationships, the ratio of 55 Co to 48 Cr can be shown, as in the 
previous section, to be proportional to the ratio of 55 Co to 52 Fe, 
with 



X( 5S Co) 
X( 48 Cr) 



1.43 XIO" 11 ^ exp (92,15) Ya 



X( 52 Fe) 



25 



X( 55 Co) 
X( 52 Fe) 



(21) 



Comparison of the right hand panel of Figs. [TT]to[9]confirms this 
approximate relationship. 
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Fig. 11. (Color online) Averaged final mass fraction ratios of vanadium to titanium and of manganese to titanium as tracers of the 
initial neutron excess. The meaning of the point types and colors, as well as the ranges of parameters are the same as in Fig. [9] Left: 
Mass ratio V/Ti. The solid line is a fit given by, M(V)/M(Ti) = 14. 8/7° 73 . The dashed line shows the relationship expected from an 
analytic model, Eq.[l9] Right: Mass ratio Mn/Ti. The solid line is a fit given by, M(Mn)/M(Ti) = 1 I7O77 068 . The dashed line shows 
the relationship expected from an analytic model, Eq. [21] 



5. Tracers of the progenitor metallicity: Optical 
epoch 

As seen in Table [2] there are several important contributors to 
the abundances of the elements from scandium to iron in the 
optical phase, most of them unstable, which implies that their 
abundances vary in time. The variation in time of their ratios 
makes it more difficult to define in this phase a tracer of the 
explosion properties. However, there are interesting correlations 
between the progenitor metallicity and several abundance ratios 
in the optical phase, namely manganese to chromium, vanadium 
to manganese, and titanium to manganese. All these ratios share 
the need to measure the abundance of manganese. We have not 
been able to found any tracer of the progenitor metallicity in this 
phase that does not involve manganese. 

Figure [12] shows the evolution of several abundance ratios 
for the extreme values of the initial neutron excess explored in 
the present work. The separation between the curves belonging 
to the different 77 makes these abundance ratios promising tracers 
of the initial neutron excess. 



5.1. Ratio of manganese to chromium 

Among the element couples shown in Fig. [12] the ratio Mn/Cr 
exhibits the best behavior, because it remains nearly constant af- 
ter day three. At these times, manganese is composed of 53 Mn, 
made as 53 Fe in the nucleosynthetic epoch, in a proportion 
higher than 80%, and the dominant isotope of chromium is 52 Cr, 
made as 52 Fe, which represents more than 90% of its abundance. 

As can be seen in Fig. [13] the ratio Mn/Cr in the optical phase 
is a good tracer of the initial neutron excess. The dispersion re- 
lated to the different parameters used in the present calculations 



of Si-b is small, although it is larger for the lowest initial neutron 
excess. The best power-law fit is 

M(Mn)/M(Cr) = 2.2rj 0A& . (22) 
5.2. Ratio of vanadium to manganese 

At the optical epoch, there are various vanadium isotopes con- 
tributing significantly to the element abundance, namely 48 V, 
49 V, and 51 V. Isotope 48 V, whose parent decays in less than one 
day, disintegrates itself with a lifetime slightly less than the time 
it takes for a typical SNIa to reach maximal brightness. Isotope 
49 V, formed from 49 Cr practically instantaneously after the ex- 
plosion, has a lifetime close to one year, which makes its contri- 
bution close to constant in the whole optical epoch. Isotope 51 V 
also forms practically instantaneously after the explosion from 
51 Mn, and decays with a lifetime of nearly one month. At the 
lowest neutron excess we considered, only 4S V contributes to 
the vanadium abundance. 

The abundance ratio of vanadium to manganese is another 
good tracer of the initial neutron excess, although it presents the 
difficulty that the ratio changes with time. The left hand panel of 
Fig.[T4lshows this ratio and a power-law fit for two times. At 15 
days, close to maximal brightness of a typical SNIa, it is given 
by 

M(V)/M(Mn) = 0.018/77 042 , (23) 

whereas at 30 days, a time at which the photosphere is close to 
the center of the ejecta, the fit is 

M(V)/M(Mn) = O.O26/77 29 . (24) 
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Fig. 14. (Color online) Averaged final mass fraction ratios V/Mn and Ti/Mn in the optical phase, as a function of the initial neutron 
excess. The point types and colors in these plots refer to the time at which the ratios are drawn, blue solid circles are calculated 
at day 15, whereas red solid squares belong to day 30. The ranges of the parameters r r j se , t, and X( l2 C) are as in Fig. [9] Left: 
Abundance ratio of vanadium to manganese. The ratios are fit by power-law functions, namely M(V)/M(Mn) = O.OI8/77 042 for the 
points at day 15, and M(V)/M(Mn) = 0.026/jf 29 for the points at day 30. Right: Abundance ratio of titanium to manganese. In 
this case, the points belonging to days 15 and 30 can be reasonably fit by a single function, although it cannot be a simple power 
law. The fit function shown is, M(Ti)/M(Mn) = exp(-6.45 - 1.1 In 77 - 0.043 In 2 77). 



5.3. Ratio of titanium to manganese 

The dominant isotope of titanium after day 10 is 48 Ti, made as 
48 Cr, whose abundance amounts to more than 90% of the total 
mass of titanium. Although the lifetime of the radioactive parent 
of 48 Ti, i.e. 48 V, is 16 days, the ratio Ti/Mn changes little after 
day ~ 15 in comparison with its dependence on the initial neu- 
tron excess. The right hand panel of Fig. [14] shows that Ti/Mn 
can be fit by the same function of 77 at 15 and 30 days, 

M(Ti)/M(Mn) = exp (-6.45 - 1 . 1 In 77 - 0.043 In 2 77) . (25) 

6. Tracers of the explosion conditions 

As seen in the previous section, the ratios of abundances of most 
iron QSE-group elements are especially sensitive to the neutron 
excess. Consequently, it is difficult to find tracers of the rest of 
parameters, namely expansion timescale, temperature rise time, 
and initial carbon abundance. We have only found one such pos- 
sible tracer, V/Mn in the X-ray epoch. 

6.1. Expansion timescale: Ratio of vanadium to manganese 
at the X-ray epoch 

The most promising pairs of elements to look for a tracer of the 
expansion timescale, t, are those that come from parents linked 
by an alpha capture in QSE, because this is the reaction type that 
depends the less in the n/p ratio, hence in the neutron excess. The 
first possibility is the pair formed by titanium and chromium, 
at the X-ray epoch, whose grandparents in the nucleosynthetic 
epoch are 48 Cr and 52 Fe, respectively. However, the results of 



the present calculations show that the ratio of their abundances 
is insensitive to r. 

The second possibility is the pair formed by vanadium and 
manganese, also at the X-ray epoch, whose grandparents are 
51 Mn and 55 Co, respectively. Their abundance ratio is only 
mildly sensitive to r and strongly dependent on 77. An approxi- 
mate measure of the expansion timescale would only be possible 
once the neutron excess is known, and only for low-metallicity 
progenitors, Z ^ O.1Z0. At higher values of the initial neutron 
excess, the ratio of vanadium to manganese abundances, V/Mn, 
is insensitive to r. 

Figure[T5]shows the ratio V/Mn for a low-metallicity progen- 
itor as a function of the expansion timescale. The calculated ra- 
tio depends on r only for expansion timescales less than ~ 0.8 s. 
The dispersion introduced by the different peak-temperature dis- 
tribution functions used is comparable to the difference in V/Mn 
for close t, so this tracer would not even allow determining the 
expansion timescale with high precision. The function, 

M(V)/M(Mn) = log [log (1 /t) 7 ] - 1 .7 , (26) 
provides a reasonable fit to the present results. 

6.2. Deflagration or detonation? 

Knowing the temperature rise timescale would allow discerning 
whether burning in the incomplete Si-b regime in SNIa proceeds 
through a deflagration or a detonation wave. Unfortunately, as 
can be expected from the discussion in section U. 1.41 we did not 
found any pair of iron QSE-group elements whose mass ratio is 
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Fig. 12. (Color online) Evolution with time of the abundance ra- 
tios of potential tracers of the progenitor metallicity in the opti- 
cal epoch: solid black lines belong to Mn/Cr, dashed blue lines 
to V/Mn, and dotted red lines to Ti/Mn. For each couple of el- 
ements, the labeled curve is the one belonging to X( 22 Ne) = 
2.5 x 10 4 , whereas the other curve with the same type and color 
belongs to X( 22 Ne) = 0.075. The dot-dashed green line shows 
the evolution of Ti/Cr ratio, which could be used as a control 
variable, see section|7] 



Fig. 13. (Color online) Averaged final mass fraction ratio Mn/Cr 
in the optical phase, as a function of the initial neutron excess. 
The meaning of the point types and colors is the same as in 
Fig. [9] The ranges of the parameters r^g, r, and X( l2 C) are as in 
Fig.0 The solid line is a fit given by, M(Mn)/M(Cr) = 2.2r] 0M . 



from these elements comes from a region processed to incom- 
plete Si-b, NSE, or a mixture of both. 



sensitive to the rise time. The same applies to the determination 
of the initial carbon mass fraction. 

7. The ratio of titanium to chromium as a control 
variable 

We have already explained that the ratio of the abundance of tita- 
nium to that of chromium, Ti/Cr, is insensitive to the parameters 
of the incomplete Si-b calculation, in the limits explored in the 
present work; i.e., it cannot be used as a tracer of any one of these 
parameters, including the neutron excess. However, this does not 
mean that measuring Ti/Cr is useless. On the contrary, it can be 
utilized as a control variable to check that the SNIa properties 
match the predictions of the present Si-b calculations. 

Figure Q~2] shows the evolution of Ti/Cr in time until two 
months after a SNIa explosion. Titanium abundance grows 
steadily as a result of the decay of 48 V, whereas the abundance of 
chromium remains stable after a few hours from the thermonu- 
clear runaway. However, the ratio Ti/Cr does not change signifi- 
cantly after day 30, and what is most important, it is practically 
the same irrespective of the metallicity of the progenitor WD. 
In the X-ray epoch, the full set of parameters used in the Si-b 
calculations leads to values of Ti/Cr in the range 0.035 - 0.06. 

The profile of Ti/Cr in the ejecta of the reference SNIa model 
of the present work, DDTc, is flat in the region undergoing in- 
complete Si-b (4.3 < Tg pgak < 5.2) with a value in the range 
Ti/Cr =i 0.025 - 0.05, whereas it rises in the inner NSE region 
by nearly an order of magnitude. As a result, measuring Ti/Cr 
might provide a way to evaluate whether the observed signal 



8. Discussion 

In this work, we investigate the properties of the nucleosynthesis 
in the incomplete Si-b regime, using a parameterized description 
of the initial chemical composition and the evolution of the ther- 
modynamic variables, in the context of SNIa. We focus the study 
on potential nucleosynthetic tracers of the parameters of Si-b. 
We discarded tracers based on the measurement of the abun- 
dance of a single element, because it would depend on the ex- 
tent of combustion within the exploding WD, which is expected 
to vary substantially from event to event. We therefore consider 
abundance ratios as the most promising targets for the present 
study. We find that the ratios that involve elements belonging to 
the silicon QSE-group do not correlate well with the parame- 
ters of Si-b, so the selected ratios are those involving the light- 
est elements of the iron QSE-group, i.e., titanium, vanadium, 
chromium, and manganese. We hereby summarize the findings 
of this work and briefly discuss the plausibility of measuring 
the most interesting abundance ratios either in SNIa remnants 
or in SNIa near maximum brightness, the so-called photospheric 
epoch. 

8.1. Young supernova remnants 

A young supernova remnant is one in which the reverse shock, 
which is responsible for heating and decelerating the ejecta, has 
not had time to reach the center of the ejecta. This time coin- 
cides more or less with the moment when the shocked mass of 
the interstellar medium equals the supernova ejecta mass. Since 
the X-ray emission of SNIa ejecta peaks at or near the reverse 
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Fig. 15. (Color online) Averaged final mass fraction ratio of 
vanadium to manganese as tracer of the expansion timescale, 
at fixed neutron excess, r\ = 2.27 x 1CT 5 . The meaning of 
the point types and colors is the same as in Fig. [9] The re- 
sults are shown for different T r j se in the range l(T 9 -0.2 s, and 
X( l2 C) in the range 0.3-0.7. The solid line is a fit given by 
M(V)/M(Mn) = log [log(l/T) 7 ] - 1.7. 

shock, it can reveal the chemical composition at different loca- 
tions within the remnant. Assuming chemical stratification of the 
ejecta, in the first thousand years or so after the explosion, where 
the precise time depends on a number of environment and explo- 
sion variables, the emission comes from material whose compo- 
sition was set in the explosion by the incomplete Si-b process. 

Measurement of the abundances of the four aforementioned 
elements in a young supernova remnant would provide different 
estimates of the initial neutron excess, related to the progenitor 
metallicity, which would allow a cross-check of the results. In 
restricted cases, it would allow roughly estimating the expansion 
timescale of the explosion: 

- Tracers of the initial neutron excess. 

1 . Mn/Cr is an excellent tracer of 77, whose best fit is given 
byEq.l 

2. V/Ti is a good tracer of 77, although with high dispersion 
at the lowest values of the neutron excess explored. The 
best fit is given by Eq.fTSI 

3. Mn/Ti is an excellent tracer of 77, whose best fit is given 
byEq.|20] 

- Tracer of the expansion timescale. 

1. V/Mn can be used to obtain a rough estimate of the ex- 
pansion timescale only if the neutron excess has been 
measured previously and if the corresponding progeni- 
tor metallicity is lower than ~ O.1Z0. In this case, the 
best fit is given by Eq.[26] 

- Control variable. 

1. Ti/Cr barely depends on any parameter of Si-b, so it can 
be used to check that the present predictions are correct 
or to check that the reverse shock is inside the incomplete 



Si-b zone. The value of Ti/Cr we found is in the range 
0.035 - 0.06. 

The mass of each one of these elements in the ejecta con- 
ditions their detectability. Although the absolute mass may vary 
from event to event, a typical mass of chromium is ~ 5 - 7 x 
10 3 M©, and the mass of the other three elements can be de- 
duced from Figs. l9land [TTH15l Titanium mass is one to two or- 
ders of magnitude less than chromium mass, whereas the mass of 
manganese oscillates between two orders of magnitude less and 
the same order as the chromium mass, depending on the neutron 
excess. Vanadium is the least prolific element, with masses two 
to four orders of magnitude less than the chromium mass. 

Chromium and manganese have alre ady been detected with 
Suzaku in a number of SNIa remnants (lYamaguchi & Kovamal 
120101 : lYamaguchi et al.ll2.Q12b . The detection of vanadium in X- 
rays is problematic due to its low abundance, whereas that of 
titanium may be made difficult by superposition of the emission 
from non-dominant ions of lavishly produced calcium, i.e. those 
emitting at a higher energy than the bulk of calcium. The emis- 
sion peaks of titanium to manganese are separated well in en- 
ergy, by more than ~ 300 eV, which ensures that their features 
can be resolved by the instrument XIS onboard Suzaku, as well 
as with future missions. For instance, the IXO proposal contem- 
plated a resolution of 6 eV at 6 keV, with an effective area of 
0.65 m 2 , which is ~ 6.5 times better than the effective area of 
XIS at the same energy. Thus, one can expect that detection of 
chromium and manganese, and perhaps titanium, in supernova 
remnants will be common in the future. 

A quite different question is the quantitative interpretation of 
X-ray observations of supernova remnants, needed to obtain the 
desired abundance ratios. First, not all the atomic data needed to 
properly model X-ray emission from these elements are avail- 
able at present. Second, it is difficult to estimate the physical 
variables that influence X-ray spectra. Hopefully, technical ad- 
vances in the near future will make this interpretation more reli- 
able. 

8.2. Photospheric epoch of SNIa 

The photospheric epoch of SNIa starts shortly after the explo- 
sion, and ends approximately a month later. It shares a common 
picture with young supernova remnant, in the sense that the opti- 
cal emission is determined mainly by the position of a receding 
front, in this case the photosphere. While the photosphere re- 
mains far enough from the center of the ejecta, the optical emis- 
sion thus comes from matter whose chemical composition was 
set by the incomplete Si-b process. In this epoch, it is necessary 
to measure the abundance of manganese to obtain any estimate 
of the initial neutron excess. 

- Tracers of the initial neutron excess. 

1. Mn/Cr again provides the best measurement of 77, be- 
cause both their ratios vary little with time, and it corre- 
lates strongly with the neutron excess. The best fit, which 
is applicable after day five after the explosion, is given by 
Eq.|22] 

2. V/Mn varies with time owing to the different contribu- 
tions of radioactive isotopes of vanadium. The best fits 
are given by Eq. [23] (day 15) and Eq. [24] (from day 30 
on). 

3. Ti/Mn can be fit as a single function of 77, Eq. [25] from 
day 10 on. 

- Control variable. 
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Ti/Cr varies with time but is insensitive to any of the param- 
eters of incomplete Si-b (Fig. [12J. 

lHatano et al.l (Q999) studied the ion signatures in SNIa spec- 
tra, concluding that there are several isotopes from titanium to 
manganese that are potentially identifiable in the photospheric 
epoch, namely Ti n, V n, Cr i, Cr n, Mn i, and Mn n. However, the 
only elements detected so far are titanium and chromium, both 
in SN 2002bo and SN2004eo dStehle et al.ll2005ar. iMazzali et alJ 
2008). Continued improvement in the techniques used in the 
recording and interpretation of SNIa spectra, will provide better 
estimates of the mass of titanium, chromium, and manganese in 
future observations^ Obviously, the next SNIa exploding in the 
Milky Way will provide the opportunity to measure its chemical 
composition with unprecedented detail. 

8.3. Final remarks 

The use of the proposed abundance ratios to estimate SNIa prop- 
erties is robust, in the sense that 

1 . The profile of each one of the abundance ratios is homoge- 
neous throughout the region affected by incomplete Si-b. As 
a result, it does not depend on the precise position of the 
reverse shock (X-ray epoch) or of the photosphere (optical 
epoch). 

2. A control variable is defined, which is the ratio Ti/Cr that, if 
measured, would allow checking that the emitting region is 
within the incomplete Si-b layers. 

3. In SNIa models, titanium and chromium are synthesized 
mainly by the incomplete Si-b process. Manganese can have 
an additional contribution from NSE regions, whose rele- 
vance depends on both the central density of the progeni- 
tor WD at thermal runaway and the initial development of 
the combustion wave. Vanadium can only take contributions 
from a slightly wider range of peak temperatures than ex- 
plored here. 

It is important that the empirical correlations between abun- 
dance ratios and initial neutron excess, found as a result of the 
integration of the nuclear evolutionary equations, are supported 
by a theoretical analysis based on the properties of QSE groups 
in Si-b. This theoretical analysis allows us to understand the con- 
ditions under which one can expect the empirical correlations to 
apply. It has also been shown that the relatively low entropy of 
incomplete Si-b in SNIa is a necessary condition for the tight 
correlation between the abundance ratio Mn/Cr and the initial 
neutron excess. This correlation therefore cannot be extrapolated 
to incomplete Si-b matter in core-collapse supernovae. 
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